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1. Introduction 

The IceTop air shower array now under con- 
struction at the South Pole as the sur face compo- 

nent of the IceCube neutrino telescope (jAchterberg et al 
20061) detected an unusual near-solar-minimum 
Ground Level Enhancement (GLE) after a solar 
flare on 13 December 2006. Beginning at 0220 
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UT, the 4B class flare occurred at solar coordi- 
nates S06 W24, accompanied by strong (X3.4) 
X-ray emission and type II and IV radio bursts. 
The LASCO coronagraph on the SOHO spacecraft 
observed a halo CME launch from the Sun at ~ 
0225 UT with sp eed estimated to be ~ 1770 km/s. 
We have begun (jBieber et al.ll2007l ) a comprehen- 
sive analysis of the propagation of solar energetic 
particles in this event. However the focus of this 
Letter is the new and unique ability of IceTop to 
derive the energy spectrum of these particles in 
the multi-GeV regime from a single detector with 
a well defined viewing direction. 

When completed, IceTop will have approxi- 
mately 500 square meters of ice Cherenkov col- 
lecting area arranged in an array of 80 stations on 
a 125 m triangular grid to detect air showers from 
one PeV to one EeV. Each station consists of two, 
two meter diameter tanks filled with ice to a depth 
of 90 cm. Tanks are instrumented with two Dig- 
ital Optical Modules (DOM) operated at differ- 
ent gain settings to provide appropriate dynamic 
range to cover both large and small air showers. 
Each DOM contains a 10 inch photomultiplier and 
an advanced readout system capable of digitiz- 
ing the full waveform. For historical reasons, the 
two discriminator counting rates recorded in each 
DOM are termed SPE (Single Photo Electron), 
and MPE (Multi Photo Electron). In the present 
analysis the SPE threshold corresponds approxi- 
mately to 20 photoelectrons (PE), and the MPE 
threshold to 100 PE. 

Due to the high altitude (2835m) and the nearly 
zero geomagnetic cutoff at the South Pole, sec- 
ondary particle spectra at the detector retain a 
significant amount of information on the spectra 
of the primary particles. In a thin, ionization de- 
tector these secondary particles either would not 
interact, or would produce virtually indistinguish- 
able signals. This is not the case in the thick Ice- 
Top detector, where a traversing muon produces 
130 PE and the typical electron only 15 PE. Signal 
amplitude therefore carries information about the 
composition and spectra of the incident particles, 
albeit integrated over broad regions of the spec- 
trum. In particular, differences in counting rates 
of discriminators at different thresholds allow us 
to infer the particle spectrum incident at the top 
of the atmosphere. 
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13 December 2006 (UT) 

Fig. 1. — Average Single Photo Electron (SPE) 
discriminator counting rate in IceTop prior to and 
during the solar particle event of 13 December 
2006. Data are averaged over five minute inter- 
vals. 
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Fig. 3. — IceTop proton spectrum from the fit in 
Figure 2 (heavy blue line with one sigma error 
band). The black line is the assumed background 
cosmic ray proton spectrum and the points are 
the maximum proton fluxes from GOES spacecraft 
data. 
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Fig. 2. — Increase of individual IceTop discrim- 
inator counting rates averaged from 0320UT to 
0410UT over that during a reference interval 
(0115UT to 0255UT) plotted against counting rate 
during the reference interval. Bars indicate statis- 
tical errors. Blue line is the best fit to a power law 
(in rigidity) spectrum, AP'' . Gray lines give the 
one sigma error band on the spectrum as discussed 
in the text. Three discriminators (open symbols) 
showed anomalous behavior and were excluded as 
described in the text. The insert illustrates cor- 
relation of error estimates for the fit parameters 
(A 7). 
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December 13, 2006 

Fig. 4. — (Top) Calculated increase in a sea level 
neutron monitor based on the IceTop spectrum 
(heavy black circles) compared with the count- 
ing rate of several near sea level neutron mon- 
itors. Traces for Oulu, Apatity, Mawson, No- 
rilsk and Barentsburg are labelled while those 
for Cape Schmidt, Fort Smith, Inuvik, Nain, 
Peawanuck, Tixie Bay, and Thule are not. (Bot- 
to m) Spectral index from Ic eTop compared to that 
of IVashenvuk et all (|2007bh 
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On 13 December 2006, IceTop was returning 
useful data from 32 SPE and 15 MPE discrimi- 
nators operating in 16 tanks. Figure 1 shows the 
average SPE discriminator counting rate as a func- 
tion of time with a pronoimccd increase due to the 
arrival of solar particles. In Figure 2 we focus on 
the time interval 0320 to 0410 where the spectral 
shape (see below) is essentially constant. We plot 
the additional counting rate (due to the solar par- 
ticles) for individual discriminators as a function 
of counting rate prior to the onset of the solar par- 
ticle fluxes. IceTop was in test mode with all dis- 
criminators set at nominal, uncalibrated thresh- 
olds. Fortuitously this produced a significant gap 
between MPE and SPE count rates, as well as a 
modest range of true thresholds within each group. 
To derive an energy spectrum from these data we 
had to deal simultaneously with the unknown el- 
ement composition of the particles (primarily the 
solar proton/alpha ratio) and the lack of a detec- 
tor calibration. 

Yield functions, with units arca-solid-angle, de- 
scribe the relation between particle flux at the 
top of the atmosphere and the occurrence rate 
of a specified signal. For IceTop they depend on 
the arrival direction, rigidity (P) and mass of the 
primary nucleus and on the discriminator setting 
which determines the light level required to count 
a particle. At high latitudes such as South Pole 
defocusing in the geomagnetic field produces an 
isotropic flux at the top of the atmosphere even 
if the flux outside the magnetosphere is highly 
anisotropic. At low energy (when the probabil- 
ity of a particle or its progeny to reach the surface 
is small), the yield function is smaller than the 
physical area-solid-angle of the tank, but at high 
energy (when a shower can give rise to a signal 
even if the trajectory of the primary passes out- 
side the tank) it is larger. 

By definition, the convolution of a yield func- 
tion Spe{P) with a particle spectrum J{P) gives 
the counting rate above the corresponding pe 
threshold. The product of the yield function and 
the spectrum is termed the response function. 
Thus, at time t during the solar flare the counting 
rate of the ith discriminator, with threshold pe(z). 



E 



5p%.)(P){Jfc(P) + AJfe(P,i)}dP, 

(1) 

where J{P) is the steady state cosmic-ray spec- 
trum, A J(P, t) is the additional flux of particles 
at time t during the event, and the summation is 
over particle species. Note that by interchanging 
summation and integration the concept of a single 
response function representing a composite spec- 
trum is well defined. 

1993[ ). and mea- 



Using FLUKA (jFasso et al 



sured cosmic ray composition and spectra appro- 
priate to solar minimum, we generated galactic 
cosmic ray (GCR) response functions for several 
different thresholds, which could be interpolated 
to produce response functions for arbitrary thresh- 
olds. Lacking a calibration, we integrated the re- 
sponse functions to predict the background count- 
ing rate for each possible threshold, then assigned 
to each discriminator the response function that 
exactly predicted the observed counting rates dur- 
ing the reference interval 0115UT to 0255UT. 
Most derivation s of a solar spectrum 
(120061 )1 



Bombardieri et al 



(.e.g. 

parti- 



assume the 

cles to be all protons and construct a proton 
yield function from GCR respo nse functions (e.g. 
Lockwood fc Debrunner ( 19991 )). Composition at 
these energies has never actually been measured, 
so we make the much simpler assumption that the 
composition is galactic, constructing yield func- 
tions by dividing the interpolated GCR response 
functions by a modifi ed force field proton spec- 
trum ( see Figure 3: Caballero-Lopez fc Moraall 
(|2004[ ): IClem et all (l200j)). If solar particles are 
proton rich compared to galactic, these yield func- 
tions produce a lower limit on the solar proton 
spectrum. We estimate that for an all proton 
composition the true intensity would be ~ 1.2 
times this limit. 

With this set of yield functions we minimized 
for a power law (in rigidity) spectrum. The 
actual fit is shown in Figure 2. The minimum was 
1.5 per degree of freedom, confirming the visual 
impression that some of the deviation from the 
fit is due to inherent differences among the detec- 
tors. Results from a Monte Carlo simulation are 
presented in an insert, with the parameter pair 
resulting from each realization plotted as a point. 
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The 68.3% of the reahzations producing the lowest 
are indicated in red; we take this to define the 
one sigma error range of the parameters. Due to 
the high correlation we do not discuss the errors on 
the parameters separately. The error band shown 
is the envelope of all possibilities within the 68.3% 
contour. Three discriminators (open symbols in 
Figure 2) were excluded from the fit primarily be- 
cause their background counting rates stand so far 
outside the cluster of discriminators operating un- 
der the same nominal conditions. Correlation of 
their background count rate with barometer read- 
ing over the three day interval surrounding the 
flare event was also anomalous compared to that 
of the other discriminators. 

The derived IceTop spectrum is shown in Fig- 
ure 3 by the blue curve, with the heavier line de- 
noting the energy range that contributes substan- 
tially (tenth to ninetieth percentile) to the fit. The 
parameter ranges quoted correspond to the dashed 
rectangle in Figure 2. The overshoot of the spec- 
tral extrapolation compared to low energy proton 
fluxes from the GOES spacecraft is rather typi- 
cal, and is generally interpreted as a steepening of 
the spectrum over the intervening energy range. 
Note that IceTop is able to derive the spectrum 
of a small increase over a large background. An 
increase of this magnitude would not be statisti- 
cally resolved with a detector of a size practical 
for flight on a spacecraft or balloon. 

3. Interpretation 

Repeating this analysis on 20 minute subsets of 
the data, and using a published response function 
(jMoraal et al.lll989l ) we calculated (Figure 4) the 
expected count rate increase for a sea level neu- 
tron monitor. The error bars were determined, as 
in Figure 2, as the range of allowed values from pa- 
rameter pairs within a 68.3% contour. If, instead, 
we were to compute the standard deviation over 
all realizations in the Monte Carlo simulation, the 
error bars would be a factor ~ 0.65 smaller. For 
comparison we show one minute averages of the 
counting rates of several near sea level neutron 
monitors. Particles first arrived from the sunward 
direction, observed best by Oulu, Mawson, and 
Apatity in a very tight b eam focused by the inter- 
planetary magnetic field teieber et al. 2007 ). The 
initial rise was seen by Mawson, then the beam 



moved over to Oulu and Apatity, before briefly 
switching back toward Mawson. By the time Bar- 
entsburg saw the beam, nearly isotropic scattered 
particles were beginning to dominate the flux. The 
viewing direction of IceTop was similar to that 
of monitors that primarily observed backscattered 
particles. We consider the agreement of our calcu- 
lation with the observations from these monitors 
to be rather good, given that we arguably have no 
free parameters. 

Traditional methods for determining energy 
spectra rely on observations from pairs or groups 
of stations with different geomagnetic cutoff and 
station altitude (' Lockwood et al.ll2002l : iRyan et"al] 



2005t feonibar dieri et all l20oi) which typically 



have strongly energy dependent viewing direc- 
tions. Such an analysis has been report ed for 



the 13 December 2006 event by Va shenvuk et al 



(l2007al lbh who employ the world network of neu- 
tron monitors to achieve a range of viewing direc- 
tions and geomagnetic cutoffs under the common 
assumption that the event can be modelled by a 
function separable in energy and anisotropy. In 
the lower panel of Figure 4 we compare the spec- 
tral index derived from both analyses. Late in 
the event, when the fluxes are more is otropic, the 



agreem ent is excellent. Early on. iVashenvuk et al 



(20073) derived a significantly harder overall spec- 
trum than ours, and also reported a pronounced 
softening of the spectrum with time. In contrast, 
our analysis yielded a nearly constant spectral in- 
dex. We believe that the d i scre pancy results from 
the way 'Vas henvuk et all (|2007b) parameterized 
the anisotropy. As a result, their fit confuses 
anisotropy evolution with spectral evolution. We 
are confident that we eventually will converge on 
a common understanding when the precise spec- 
trum derived from IceTop is properly included as a 
constraint on the fit. When definitive results from 
the PAMELA spacecraft instrument are available 
they should also contribute greatly to a compre- 
hensive analysis. It is also clear that future results 
from IceTop will be greatly enhanced by the neu- 
tron monitor network, which will continue to be 
the primary source of information on anisotropy. 

4. Future Plans for IceTop 

Encouraged that such a straightforward analy- 
sis of IceTop data yields a useful picture of the 
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time dependent spectrum of the solar particle 
event on 13 December 2006, we are working to 
better understand the instrumentation to reduce 
systematic uncertainties. Wc arc also reconfigur- 
ing IceTop to increase statistical precision. DOMs 
are being reprogrammed to collect and transmit 
histograms, with ten second resolution, of the in- 
tegrated charge signal for all events that trigger 
the MPE discriminators. These will then be set to 
trigger at a rate of approximately 2000 Hz, a limit 
determined by the tolerable dead time of the sys- 
tem for air shower studies. Each of the (eventual) 
160 high gain DOMs will thus return a spectrum 
that is statistically equivalent to the spectrum re- 
turned by the entire ensemble of DOMs employing 
the present data collection method. SPE discrim- 
inators will be set at a variety of thresholds to 
populate the regime at rates higher than 2000 Hz, 
probably up to about 10,000 Hz. DOMs are ca- 
pable of accumulating histograms at these higher 
rates, but the dead time would be unacceptable 
for normal operation. Overall wc should be able 
to achieve a sensitivity two orders of magnitude 
greater than current neutron monitors for event 
detection. 

These changes will dramatically improve the 
energy resolution of IceTop for determining solar 
flare spectra. In Figure 2, there is some spread 
within the two clusters of points but not enough to 
go much beyond a two parameter fit. With prop- 
erly spaced coverage and extensions out to 10,000 
Hz, we will be able to measure spectral curvatures 
or cutoffs. We are indeed looking forward to the 
coming solar maximum. 
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